Quantum Monte Carlo study of the 3D-1D crossover for a trapped Bose gas 
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We investigate the cross-over from three to one dimension in a Bose gas confined in highly 
anisotropic traps. By using Quantum Monte-Carlo techniques, we solve the many-body Schrodinger 
equation for the ground state and obtain exact results for the energy per particle and the mean 
square radii of the cloud in the transverse and longitudinal direction. Results are compared with 
the predictions of mean-field theory obtained from the Gross-Pitaevskii equation and with the ID 
Lieb-Liniger equation of state. We explicitly prove the occurrence of important beyond mean-field 
effects, including the appearance of Fermi-like properties as the system enters the Tonks gas regime. 

PACS numbers: 03.75.Fi, 05.30.Fk, 67.40. Db 



I. INTRODUCTION 

The study of trapped Bose systems in low dimensions 
is currently attracting a lot of interest. In particular, 
ID systems are expected to exhibit remarkable proper- 
ties which are far from the mean-field description and are 
not present in 2D and 3D. The peculiarity of ID physics 
consists in the role played by fluctuations, which destroy 
long-range order even at zero temperature Q, and in 
the occurrence of characteristic effects due to correlations 
such as the fermionization of the gas in the Tonks regime 
0. Recent experiments with highly anisotropic, quasi- 
one-dimensional traps have shown first evidences of ID 
features in the aspect ratio and energy of the released 
cloud as well as in the coherence properties of con- 
densates with fluctuating phase From a theoretical 
viewpoint, the emergence of ID effects in the properties 
of binary atomic collisions, by increasing the confinement 
in the transverse direction, has been pointed out in j(| . In 
the case of harmonically trapped gases, the occurrence of 
various regimes possessing true or quasi-condensate and 
the possibility of entering the Tonks gas regime of im- 
penetrable bosons has been discussed in |7|- 

The ground-state properties and excitation spectrum 
of a homogeneous ID system of bosons interacting 
through a repulsive contact potential have been calcu- 
lated exactly by Lieb and Liniger long time ago ||. For 
a fixed interaction strength the Lieb-Liniger equation of 
state reproduces in the high density regime the mean- 
field result obtained using the Bogoliubov model and 
in the opposite limit of low density coincides with the 
ground-state of impenetrable bosons [||. For ID sys- 
tems in harmonic traps, the exact many-body ground- 
state wavefunction in the Tonks regime has been recently 
calculated |9), and the equation of state interpolating 
between the mean-field and the Tonks regime has been 
obtained within the local density approximation in JTo| . 
Methods based on local density approximation in the lon- 
gitudinal direction and on the Gross-Pitaevskii equation 
for the transverse direction have been recently employed 
to predict the frequency of the collective excitations |ll]] 



and the ground-state energy in the 3D-1D cross-over as 
well as in the ID mean field - Tonks gas cross-over Jl^] . 

In this work we present exact Quantum Monte-Carlo 
results for the 3D- ID cross-over in harmonically trapped 
Bose gases. As a function of the anisotropy parameter 
of the trap we calculate the ground-state properties of 
the system and for highly anisotropic traps we point out 
the occurrence of important beyond mean-field effects in- 
cluding the fermionization of the gas. 



II. THEORY 



We consider the following Hamiltonian 
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describing a system of N spinless bosons of mass m inter- 
acting through the two-body interatomic potential V(r) 
and subject to the axially symmetric harmonic external 
field V ex t(r) = m(o;^r^ + w^z 2 )/2, where z is the axial 
coordinate, r±_ is the radial transverse coordinate and u> z , 
lu±_ are the corresponding oscillator frequencies. 

For the interatomic potential we use two different re- 
pulsive model potentials: 
1) Hard-sphere (HS) potential defined as 



V(r) 



+oo (r < a) 
(r > a) 



(2) 



where the diameter a of the hard sphere corresponds to 

the s-wave scattering length. 

2) Soft-sphere (SS) potential defined as 



V(r) 



V 




(r < R) 
(r > R) 
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The s-wave scattering length is given by a = R[l — 
tnnh(K R)/K R], with Kg = V Q m/h 2 and V > 0. For 
finite Vq one always has R > a, while for Vq — > +oo 
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the SS potential coincides with the HS one with R = a. 
The height Vq of the potential is fixed by the value of 
the range R in units of the scattering length, for which 
we choose R = 5a. It is worth noticing that the HS 
and the SS model with R = 5a represent two extreme 
cases for a repulsive interatomic potential. In the HS 
case, the energy is entirely kinetic, while for the SS po- 
tential a ~ (m/h 2 ) J °° V(r)r 2 dr, according to Born ap- 
proximation, and the energy is almost all potential. By 
comparing the results of the two model potentials we can 
investigate to what extent the ground-state properties of 
the system depend only on the s-wave scattering length 
and not on the details of the potential. 

The relevant parameters of the problem are the num- 
ber of particles N, the ratio a/aj_ of the scattering 
length to the transverse harmonic oscillator length a± = 
yjfi/muj±_ and the anisotropy parameter A = u) z /uj±. For 
a given set of parameters we solve exactly, using the 
Diffusion Monte-Carlo (DMC) method, the many-body 
Schrodinger equation for the ground state and we calcu- 
late the energy per particle and the mean square radii 
of the cloud in the axial and radial directions. Impor- 
tance sampling is used through the trial wavefunction 
V>t(R) = ^T(ri,..,r w ) = liih( r i)Ui<jMrij), where 
fi and fi are respectively one and two-body Jastrow 
factors. For the one-body term, which accounts for the 
external confinement, we use a simple gaussian ansatz 
fi(r±,z) — exp(— ol\j\ — a z z 2 ), with a± and a z opti- 
mized variational parameters. The two-body term ft(r) 
accounts instead for the interparticle interaction and is 
chosen using the same technique employed in Ref . |l3| for 
a homogeneous system. Of course, since DMC is an exact 
method, the precise choice of ^(R) is to a large extent 
unimportant and the results obtained are not biased by 
the choice of the trial wavefunction |L4j . 

The DMC results are compared with the predictions of 
mean-field theory which are obtained from the stationary 
Gross-Pitaevskii (CP) equation 

(-^V 2 + V ext (r) + g(N - l)|$(r)| 2 ) $(r) = M $(r) , 

(4) 

where $(r) is the order parameter normalized to unity: 
J dr|<I>(r)| 2 = 1 and g = ATifi 2 a/m is the coupling con- 
stant. Further, finite size effects have been taken into 
account in the CP equation by the factor N — 1 in the 
interaction term |l5f| . In the case of anisotropic traps 
with A < 1, the GP equation (|4|) is expected to pro- 
vide a correct description of the system if the transverse 
confinement is weak a/aj_ <C 1, and if the mean sepa- 
ration distance between particles is much smaller than 
the healing length 1/n 1 / 3 <C £, where £ = \/\f%ima and 
n is the central density of the cloud. In terms of the 
linear density along z, ni£>{z) ~ 2ir J Q dr± rj_n(r±, z), 
this latter condition reads 1/riiD <C a\ja. If the mean 
separation distance between particles in the longitudinal 



direction becomes much larger than the effective ID scat- 
tering length given by a\ja 0, the mean-field approxi- 
mation breaks down because of the lack of off diagonal 
long range order. 

The system enters the ID regime when the motion in 
the radial direction becomes frozen. In this regime the ra- 
dial density profile of the cloud is fixed by the harmonic 
oscillator ground state, resulting in a mean square ra- 
di us wh ich coincides with the transverse oscillator length 
\J (r'j_) = aj_. Further, the energy per particle is dom- 
inated by the trapping potential and one has the con- 
dition E/N — huj± <C Tiuj±. If the discretization of lev- 
els in the longitudinal direction can be neglected, i.e. if 
E/N — hu!± ^S> fitu z , the ID system can be described 
within the local density approximation (LDA). In this 
case, the chemical potential of the system is calculated 
through the local equilibrium equation 

TTl 

yL = Tiw±_ + piocai{ni D {z)) + , (5) 

where Tiuj^ is the dominant contribution of the trans- 
verse confinement and piocaiinio) is the chemical po- 
tential corresponding to a homogeneous ID system of 
density n\u- If the ratio a/a± <C 1, the local chem- 
ical potential can be obtained from the Lieb-Liniger 
(LL) equation of state with the effective ID coupling 
constant g 1D = g/{2-Ka 2 A _) §,0. One finds: m oca i = 
d[n 1D e LL (n 1D )]/ dni D , where e LL is the LL energy per 
particle. By using Eq. (^) and the normalization condi- 
tion f^^dz ni£)(z) — N, one can obtain the chemical 
potential /iasa function of N . The ground-state en- 
ergy of the system with a given number of particles can 
then be calculated through direct integration of p(N). 

If n\Da\ja ^> 1, the system is weakly interacting and 
the LL equation of state coincides with the mean-field 
prediction: €ll = 5i_d"-il>/2. In this regime, one finds 
the following results for the energy per particle 

and for the mean square radius of the cloud in the longi- 
tudinal direction 

In the opposite limit, niua\/a <C 1, the system enters 
the Tonks regime and the LL equation of state has the 
Fermi-like behavior ell — it 2 h 2 n 2 D /6m. The energy per 
particle and the mean square radius of the trapped sys- 
tem are given in this case by the following expressions 



— - huj ± = —hu ± , yjjf) = J —a ± . (8) 

In terms of the parameters of the system, the two regimes 
can be identified by comparing the corresponding en- 
ergies. The mean-field energy becomes favourable if 
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N\a\/a 2 ^> 1, whereas the Tonks gas is preferred if the 
condition NXa 2 ± /a 2 <§; 1 is satisfied. 

In order to account for effects beyond local density ap- 
proximation we have also applied the DMC method to 
a system of iV particles interacting through the Lieb- 
Liniger Hamiltonian in the presence of harmonic confine- 
ment 

h 2 N d 2 

H LL = Nhuj ± - — 22 + 9id S ( Zi ~ z ^ 

i—l % i<j 

N 11 

■ (9) 

i=l 

When the number of particles is large, the properties of 
the ground state of the Hamiltonian (|9j) coincide with the 
ones obtained from the LL equation of state within LDA. 
However, for small systems one expects deviations and 
the DMC method provides us with a powerful tool. The 
relevant parameters are the same as for the 3D simula- 
tion with the Hamiltonian (^): the number of particles N, 
the ratio a/a± fixing the strength of the contact potential 
gin and the anisotropy parameter A fixing the strength 
of the longitudinal confinement in units of fiuj±. Impor- 
tance sampling is realized through the trial wavefunction 
i/j t (z 1 , ..,z N ) =H i gi(z i )ll i<j g2(z i -Zj). The one-body 
term is of gaussian form gi(z) — exp(— a z z 2 ) and the 
two-body term has been chosen as g%{z) — cos[fc(|z| — L)] 
for \z\ < L and gi{z) = 1 for \z\ > L. The condi- 
tion fctan(fcL) = a/a 2 ^ fixes the proper boundary con- 
ditions of the two-body problem in z = related to the 
(5-potential. The parameters a z and L are variational pa- 
rameters. We notice that our choice of the two-body Jas- 
trow factor reproduces both the weakly interacting and 
the Tonks regime. In fact, if kL is small, g%{z = 0) ~ 1 
and the contact potential is almost transparent. On the 
other hand, if kL approaches 7r/2, g^iz — 0) ~ and the 
contact potential behaves as an impenetrable barrier. 

III. RESULTS 

We first consider a system of very few particles (N = 5) 
and we consider different values of the ratio a/a±. Figs. 
1-2 refer to a/a± = 0.2, and we present results for the 
energy per particle and the mean square radius of the 
cloud in the longitudinal direction as a function of the 
anisotropy parameter A = u> z /uj±_ . Results from the 
GP equation ([|) and from the Lieb-Liniger equation of 
state in LDA are also shown. We find that the HS and 
SS potential give practically the same results even for 
the largest values of A, showing that for these parame- 
ters we are well within the universal regime where the 
details of the potential are irrelevant. For large val- 
ues of A the DMC results agree well with the predic- 
tions of GP equation. By decreasing A beyond mean- 
field effects become visible and both the energy per par- 
ticle and the mean square radius approach the LL result 



when NXa^/a 2 ~ 1, corresponding to A ~ 10~ 2 . Fi- 
nally, for the smallest values of the anisotropy parameter 
(A ~ 10 -4 ) we find clear evidence of the Tonks gas be- 
havior both in the energy and in the shape of the cloud. 
It is worth stressing that beyond mean-field effects oc- 
curring in the small A regime can be only obtained by 
using DMC. A Variational Monte-Carlo (VMC) calcu- 
lation based on the trial wavefunction ^y(R) described 
above, would yield results in good agreement with mean- 
field over the whole range of values of A. DMC results 
using the Lieb-Liniger Hamiltonian Hll of Eq. (^) are 
also shown and coincide with the results of the 3D Hamil- 
tonian (0) . This shows that the 3D interatomic potential 
is correctly described by the ID 5-potential even for the 
largest values of A. In fact, due to the small number of 
particles, the density profile of the cloud in the transverse 
direction is correctly described by the harmonic oscillator 
ground-state wavefunction (see Fig. 9). The ID charac- 
ter of the system is also evident from Fig. 1 which shows 
that E/N — fiuj± is always smaller than the transverse 
confining energy. Deviations of DMC results from the LL 
equation of state arise because of finite size effects. These 
effects become less and less important as A decreases and 
one enters the regime (E/N — hcu±)/frui± 3> A where 
LDA applies. In terms of the mean square radius of the 
cloud (see Fig. 2), the condition of applicability of LDA 
requires (z 2 ) 1 / 2 much larger than the corresponding ideal 
gas (IG) value. 
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FIG. 1. Energy per particle as a function of A. DMC re- 
sults: HS potential (solid circles), SS potential (solid trian- 
gles), LL Hamiltonian Q (open circles). Dashed line: GP 
equation (^), solid line: LL equation of states in LDA, dot- 
ted line: gas of Tonks, dot-dashed line: non-interacting gas. 
Error bars are smaller than the size of the symbols. 

In Figs. 3-4 we present results for a/a± = 0.04, corre- 
sponding to a less tight transverse confinement or, equiv- 
alently, to a smaller scattering length. By decreasing 
a/aj_ we enter more deeply in the universal regime where 
the theory of pseudo-potentials applies and we find no 
difference between HS and SS results. Further, as for 
the a/a± = 0.2 case, we find no difference between DMC 
results obtained starting from Hamiltonian (pi) and from 
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the ID Hamiltonian (^). Compared to Figs. 1-2, the 
cross-over between mean field and Lieb-Liniger occurs 
for smaller values of A . For the smallest values of A 
beyond mean-field effects become evident, though one 
would need to decrease A even further to enter the Tonks 
gas regime. 
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FIG. 2. Mean square radius along z as a function of A. 
Symbols and line codes are the same as in Fig. 1. Error bars 
are smaller than the size of the symbols. 
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FIG. 3. Energy per particle as a function of A. DMC re- 
sults: HS potential (solid circles), LL Hamiltonian (^) (open 
circles). Line codes are the same as in Fig. 1. 
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FIG. 4. Mean square radius along z as a function of A. 
Symbols and line codes are the same as in Fig. 3. 

The results for a/a± = 1 are shown in Figs. 5-6. In 
this case the HS and SS potential give significantly dif- 
ferent results in the large A regime. For both potentials 
the mean-field description is inadequate. By decreasing A 
the HS system enters the Tonks regime before approach- 
ing the LL results, whereas the SS system crosses from 
the ideal gas (IG) regime to the Tonks regime. For this 
value of a/a± the DMC results with the Hamiltonian (|^) 
coincide exactly with the ones of the LL equation of state 
in LDA due to the large coupling constant. 
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FIG. 5. Energy per particle as a function of A. Symbols 
and line codes are the same as in Fig. 1. 
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FIG. 6. Mean square radius along z as a function of A. 
Symbols and line codes are the same as in Fig. 1. 

Figs. 7-8 refer to a much larger system with N = 100 
and a/a± = 0.2. In this case, we see a clear cross-over 
from 3D mean field, at large A, to ID LL at small A. 
Important beyond mean-field effects become evident in 
the energy per particle as N\a\/a 2 ~ 1, corresponding 
to A ~ 10~ 3 . The Tonks regime would correspond to 
even smaller values of A which are difficult to obtain in 
our simulation. However, for the smallest values of A 
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reported in Fig. 7 we find already very good agreement 
with the LL equation of state. One should notice that 
small deviations from mean field are also visible for A ~ 
1, and are due to high density corrections to the GP 
equation. The DMC results with the ID Hamiltonian 
(Q) follow exactly the LDA prediction showing that the 
deviations seen in Figs. 1-2 are due to finite size effects. 
In the cross-over region from the mean- field to the ID 
LL regime, residual 3D effects are still present (see Fig. 
9) and produce small deviations from the LL equation of 
state. 
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FIG. 7. Energy per particle as a function of A. Symbols 
and line codes are the same as in Fig. 1. 



find */ — a ± over the whole range of values of A. 
It is worth noticing that for the N — 5 system with SS 
potential, the largest deviations from VvT) = a± are 
achieved for a/a± = 0.2, corresponding to a transverse 
confinement a± = R where R is the range of the SS po- 
tential. 
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FIG. 9. Mean square radius in the radial direction as a 
function of A. Solid symbols: HS potential; open symbols: SS 
potential. Down triangles: N = 100 and a/a± = 0.2; circles: 
N = 5 and a/a± — 1; squares: N — 5 and a/a±_ = 0.2; up 
triangles: N = 5 and a/a± = 0.04. Error bars are smaller 
than the size of the symbols. 



IV. CONCLUSIONS 
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FIG. 8. Mean square radius along z as a function of A. 
Symbols and line codes are the same as in Fig. 1. 

Finally, in Fig. 9, we show results for the mean square 
radius in the transverse direction. The cross-over from 
3D to ID is clearly visible in the case of TV = 100, for 
both the HS and SS potential, and for the HS potential 
in the case of N = 5 and a/a± = 1. For the system 
with N — 5 and a/a± = 0.2 we only see small devia- 
tions from = a± for the largest values of A. In 
the a/a± = 0.04, as well as in the a/aj_ = 1 case with 
the SS potential, the transverse density profile is well de- 
scribed by the harmonic oscillator wavefunction and we 



In this paper we present exact Quantum Monte Carlo 
results of the ground-state energy and structure of a Bose 
gas confined in highly anisotropic harmonic traps. Start- 
ing from a 3D Hamiltonian, where interparticle interac- 
tions are model by a hard-sphere or a soft-sphere poten- 
tial, we show that the system exhibits striking features 
due to particle correlations. By reducing the anisotropy 
parameter A, while the number of particles N and the ra- 
tio a/a± of scattering to transverse oscillator length are 
kept fixed, the system crosses from a regime where mean- 
field theory applies to a regime which is well described 
by the ID Lieb-Liniger equation of state in local density 
approximation. In the cross-over region both theories fail 
and one must resort to exact methods to account prop- 
erly for both finite size effects and residual 3D effects. 
For very small values of A we find clear evidence, both 
in the energy per particle and in the longitudinal size 
of the cloud, of the fermionization of the system in the 
Tonks-Girardeau regime. 
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Note added. While this work was being prepared for 
publication, a preprint by D. Blumc fig] appeared in 
which the author also calculates the ground-state prop- 
erties of Bose gases in highly anisotropic traps using Dif- 
fusion Monte Carlo techniques. The results obtained are 
similar to ours. 
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